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Abstract. In a quantum measurement, a coupling g between the system S and the apparatus A triggers the establishment of 
correlations, which provide statistical information about S. Robust registration requires A to be macroscopic, and a dynamical 
symmetry breaking of A governed by S allows the absence of any bias. Phase transitions are thus a paradigm for quantum 
measurement apparatuses, with the order parameter as pointer variable. The coupling g behaves as the source of symmetry 
breaking. The exact solution of a model where S is a single spin and A a magnetic dot (consisting of A' interacting spins and 
a phonon thermal bath) exhibits the reduction of the state as a relaxation process of the off-diagonal elements of S + A, rapid 
due to the large size of N. The registration of the diagonal elements involves a slower relaxation from the initial paramagnetic 
state of A to either one of its ferromagnetic states. If g is too weak, the measurement fails due to a "Buridan's ass" effect. The 
probability distribution for the magnetization then develops not one but two narrow peaks at the ferromagnetic values. During 
its evolution it goes through wide shapes extending between these values. 



QUANTUM MEASUREMENTS, A PROBLEM OF STATISTICAL MECHANICS 

A quantum measurement has well known features which relate the initial state of the tested system S and the final state 
of the compound system S + A including S and the apparatus A [1]. Von Neumann's reduction of the state postulates 
that an ideal measurement erases the off-diagonal blocks of the density matrix of S + A in a basis where the measured 
observable is diagonal. Born's rule provides the probability of observing on the apparatus such or such value, and this 
value is correlated with the state of S after measurement. However, in order to understand how these features arise in 
an actual experiment, we need to analyze the dynamical process undergone by the coupled system S + A, to dig out 
the time scales involved, and to exhibit the specific properties of the quantum system A required so that it can be used 
as an apparatus. 

At the initial time of the measurement, S and A are uncorrelated, with density operators r (0) and ^ (0), respectively. 
The full density operator S{t), initially equal to S{0) — r (0)^^(0), evolves according to the Liouville-von 
Neumann equation. We wish to explain how it eventually reaches at the final time ff the form 

^(ff) = ^(n,f(o)n,)'^^,- (1) 

which embeds the standard properties of ideal measurements recalled above. We denote by H,- the projection operator 
over the eigenspace of the measured observable A labelled by its eigenvalue A,-, and by the set of possible final 
states of the apparatus A. Each one, is characterized by some value of a pointer variable of A, in one-to-one 
correspondence with the value A, for S, and which can be observed or registered. The expression Q means that, in 
a set of repeated experiments, a well defined outcome / occurs simultaneously for A (on which it can be observed) 
and for S, which is thereby prepared in the state n,r (0) fl,/ Trs n,r (0) through selection of the outcome / for A. Such 
events occur with a classical probability Trs n,r (0). 

The quantum system S is microscopic, and we wish its interaction with A to be sufficiently weak so that the final 
state involves no other correlation between S and A than the ones exhibited in eq. Q. The system A should therefore 
be able to switch from ^ (0) to any one of the states under the effect of a very small perturbation. It is thus subject 
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to a bifurcation. The state (0) of A should be metastable when the interaction with S is not yet turned on, and 
it should evolve towards any one of the states when triggered by this small interaction. Moreover, a robust and 
permanent registration excludes the possibihty of transition from one state to another over a reasonable delay. Such 
states should thus be stable, and should have no overlap with one another, which imposes the apparatus A to be a 
macroscopic object. There should also exist some symmetry between the possible pointer variables, so that no bias is 
introduced by the apparatus. Otherwise the probability of observing the result i would not be Trs n,r (0) as stated by 
Bom's rule, but would also depend on the lack of symmetry of A. 

A natural means for achieving these requirements is to choose A as a macroscopic system which can undergo a 
phase transition with broken invariance. In some way or another, many actual apparatuses rely on such a transition 
(although bifurcations without broken synmietry are currently used). The pointer variable is then identified with the 
order parameter, a quantity that can be observed without perturbing the apparatus. The initial state (0) is chosen 
as the symmetric state which is metastable upon sudden cooling, and the states Mi are the possible stable states with 
broken synmietry. Whereas in statistical mechanics the invariance is broken either spontaneously or under the effect 
of a scalar source, in a quantum measurement the breaking is triggered by coupling with S. Although this coupling 
should be weak, its effects on A are considerably amplified, due to the metastability of ^(0). On the other hand, we 
wish the measurement to perturb the microscopic system S as little as possible in spite of the macroscopic nature of 
A: we wish the diagonal blocks tlif (0) n, to be kept unchanged, while the laws of quantum mechanics should prevent 
the off-diagonal blocks of f (0) to subsist. Thus, although S behaves as a source for the phase transition, its quantum 
nature is essential for a full description of the dynamics. 

From the large size of A and the small size of S we can anticipate that this dynamics of the coupled system S + A 
during the measurement involves several time scales. The variables associated with S are expected to vary rapidly from 
the outset; during this stage A hardly changes. The establishment of correlations and the registration of the result in 
the form of a finite order parameter will require a much longer delay, which is required to let the macroscopic system 
A reach equilibrium. 



A SOLVABLE MODEL 

We have worked out the above ideas on an exactly solvable model, where S is a single spin j represented by the Pauli 
operators Sx, Sy, s^, and where A simulates a magnetic dot containing spins a ^"\n = 1, 2, ... A^) which can undergo 
an Ising phase transition [2]. (For other models, see [3-8].) The order parameter is the average magnetization along z, 

N 



'?i=T7l«^n (2) 
n=l 

which takes at equilibrium the paramagnetic value (m) = above the Curie temperature, or either one of the two 
ferromagnetic values (m) = imp below. In both cases the statistical fluctuation is smaU as 1 / ^/N. Thermal equilibrium 
of A = M -I- B is ensured by a weak interaction of these magnetic degrees of freedom M with a thermal bath B, 
represented by a model describing the phonons in the dot. The fuU Hamiltonian of S -|- A has therefore the form 
H = HsA + Ha, where 

N 

HsA = -gsz £ = -NgSzm (3) 

n=l 

accounts for the interaction between the z-component of the spin s and the apparatus, and where = + ^mb + 
includes the Ising interaction 

«,«' 

between the spins of the apparatus (which all interact pairwise, unlike in our previous studies, where quartets interact), 
the couphng 

^MB = x/rI I (5) 

n=\a=x,y,z 

of these spins with the bath through phonon operators 5^ and the phonon Hamiltonian H^. The initial state 
^ (0) = (0) ^B (0) of A factorizes into contributions of the magnet and of the bath, where (0) = ^/2^ 



represents a completely disordered paramagnetic state, prepared by bringing M at a high temperature before the 
measurement, and where (0) °= e^^s/r describes equilibrium of the phonons at the temperature T . This temperature 
lies below the Curie temperature {T < J) so that ,^(0) is metastable: it may transit towards the stable ferromagnetic 
states, in a rather large time if the magnet-bath coupling 7 is weak. The bath Hamiltonian (t^ will enter the problem 

only through the autocorrelation function ^fii"' {t)B^^ ^ (f')^' which under rather general conditions has the form 

Tr^B (0)£i"^ (0^1"'^ {t') = dnn'daa'K [t - t') , (6) 

K[a^)^j_^dte-'K{t) = —^,^^. (7) 

We have denoted the Debye cutoff as F. 

Since 7 is weak, the bath can be eliminated by means of a second-order treatment of i^MB- From the Liouville-von 
Neumann equation for f^, we thus find, for the partial trace D = Ti^ ^ which describes the joint evolution of S and M 
(in the presence of the bath), the equation 



^ = ^[HsA + HM,b]+-^j: L dt' a^"\Dit)K{-t')d!,"'{t')-K{t')d^r{t')Dit) 



(8) 



We have introduced, in the interaction representation, 

d^"\t)^U{t)d^,"^uHt) , (9) 

U{t)^exp[{HsA+HM)t/ih] . (10) 

DISAPPEARANCE OF SCHRODINGER CATS 

During an ideal quantum measurement, the quantity to be measured should not change. This is expressed here by the 
commutation of H with or equivalently with the projection operators FI, on the eigenstates / =1 or | of Sj. Hence, eq. 
(|8j can be decomposed into four equations governing the blocks Dij = TliDYlj —1^ or J,). Each Dij is an operator 
in the space of M. 

This decoupling allows us to treat separately the off-diagonal and diagonal blocks. The evolution of — 
governed by (|8} has been studied elsewhere [2, 9]. We briefly recall the results here. During a very brief reduction 
time Tred = fi/V^Ng, all elements of which describe correlations between the components Sx or Sy of the spin S 
and the various spins ff'"' of A, decrease down to zero. The initial order exhibited by the non-vanishing value of (s^) 
or (^Sy') is scattered into a very large number of small correlations between Sx or Sy and the z-component of the many 
spins of A. 

This relaxation process is governed by the interaction term H^j^. If the Hamiltonian if S + A did reduce to this 
term, the lost order would surge back in at the time nTi/llg), producing detectable effects through a recurrence 
process. Such recurrences are hindered, and irreversibility of the collapse is ensured, owing to the presence of the 
spin-bath coupling term Hmb. provided 7 ^ /Nh^T^. The initial order is then lost once and for all into the bath 
degrees of freedom. In the language of NMR, this irretrievable loss of phase coherence is similar to a spin-lattice 
relaxation process. However, an alternative mechanism involving no phonon bath can also suppress the recurrences, 
for an apparatus A consisting only of the magnet M. Indeed, all the matrix elements of D|j remain negligible at 
all times after Tred if the system-magnet interaction i/sA involves slightly different coupling constants gn between 
and the apparatus spins di"^, such that their relative fluctuation satisfies 5g/g 3> This second mechanism is 

comparable to the relaxation in NMR due to inhomogeneity of the external Larmor field. As in Hahn's spin echoes, 
one can imagine to retrieve here the initial order associated with Z)|| by means of an adequate setup of pulses. 

Altogether, the von Neumann collapse of the state, which can be regarded as the disappearance of "Schrodingercats" 
in a measurement, is explained as a relaxation process, which is rapid due to the large size of A^. An important role is 
also played by the mixed nature of the initial state of our apparatus. In most theoretical discussions of the measurement 



process a pure initial state is assumed. This is not realistic, however, since it is impossible to control the macroscopic 
number of degrees of freedom of an apparatus. The relaxation related to the collapse should not be confused with 
standard decoherence. It does not necessarily require the thermal bath, it depends on the measured observable, and it 
takes place over a time Vd which involves g, not T. 



REGISTRATION BY THE APPARATUS 



We now focus on the evolution of the diagonal blocks and D^, which describe the expectation values of and of 
the spin operators of M, their fluctuations and all their correlations. This evolution should end up in a final state of the 
form Q, expressing the registration of the measurement: a complete correlation established between the final value ±1 
of and the final ferromagnetic equilibrium state, characterized by the sign of the order parameter (m) = imp. The 
diagonal block Z)|| is associated with the occurrence of = +1. We have to prove that in this block the probability 
distribution for m tends for large times to become sharply peaked around +mp, without any contribution from the 
region of —mp. In the equation of motion for obtained from (|8|i, the operator i/sA depends only on the 

eigenvalue of s, in this sector, equal here to +1, and on the observable m, the successive eigenvalues of which are 
separated by a distance 2/N. This introduces through 0, and jlO> the operators 
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= ^2{g + Jm) - 2J/N = T2h [m) - 2J/N ., 



where h (in) = g+ Jm behaves as an operator-valued self-consistent field. We also introduce the function 

sm{(o' -(o)t 



K,{co)= / dse-'""K{s) = 



dco'Kico' 



n{(o' -co) 

in which (o can be replaced by the operators Ci±, and thus find the reduced equation of motion 
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(11) 



(12) 



(13) 



This evolution conserves the trace TrM^n = Trsn|r(0) = '"ii(O), a normalization consistent with Born's rule. 
Moreover the operator in the space of M turns out to be simply a function of m at each time, (f ) = A|| (m, f ), 
since this property holds at the initial time and is preserved by the motion J13t . In fact, the knowledge of A|| (m,f) 
is equivalent to that of the conditional probability {m,t) for m to take the discrete values m = — 1, — 1 +2/N, 
I —2/N, 1 if ij; equals +1, which is expressed by 



Pd{m,t) 



N\ 



A|| (m,f) 



[iA^(l+m)]![iA^(l 
From il3\ we find the equation of motion for {m, t), 
dPd{m,t) jN 



! TrD 



TT 



dt 
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N' 



-Kr{a+){l-m)Pa{m,t) 



(14) 



(15) 



-Kr ( 1 - m + Pd (m - ] - K, (£2_) (1 +m)Pd M 



where fi± are functions of m defined by ea. Jl 11 1. For shorthand we have denoted by fd {iii,t) instead of P| {m,t) the 
conditional probability of m associated with the value i- = +1 of the measured observable of S, while the subscript 
'd' indicates that m is discrete here. We can introduce likewise for the sector |J, the conditional probability P^ {m,t) 
associated with = — 1, which is obtained by merely changing g into —g in the expression il 1> of D,±. 

The registration process of the measurement is therefore just the same, in the sector t|, as the relaxation of the Ising 
model towards equilibrium under the influence of a field +g and of the phonon bath, a problem that we now study. 
The dynamical equation ( I15> has been derived from the Hamiltonian //a ^ Ngm without any other approximation than 



a weak magnet-bath coupling 7^1. The bath occurs through Kt defined by Q, ( I12t . while the Hamiltonian of M 
(including the field +g) occurs through the energy shifts hQ.± given by (^Ji- 

For sufficiently small 7 the evolution is slow and its time scale is large compared to Ti/T, so that Kt (c) can 
be replaced in il5i by K{a)). In such a short-memory approximation, we have K{—a)) = K{(£))e^^'^l^ , and the 
equation M5\ for {m,t) can be identified with a balance equation, where the probability of each spin flip induced 
by the coupling with the bath is given by Fermi's golden rule, and which might have been written directly on 
phenomenological grounds. In terms of the entropy S and the energy t/ of M associated with the density operator 
eq. il5\ satisfies in this regime an //-theorem 
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(1 +m)Pa{m) 



(16) 



which implies that Djj tends within normalization to the equilibrium distribution expA^ (^gm + jJm^) /T, with possible 
invariance breaking for small g. 



DYNAMICS OF THE PHASE TRANSITION 



Our purpose is to study the dynamics of this relaxation, which for T < J involves a bifurcation from the vicinity of 
m — towards that of imp. In the large limit, we may treat m as a continuous variable, P{m,t) — {N /2)Pa {m,t) 
being now normalized as J^^ dmP(jn,t) = 1. From (0, il l\ and il5\ we get for t ^ h/T, keeping the terms of order 
1 and 1 /N, 

dP{m,t) d . ^p. ,^ 
— 'V[m)P(m,t)- 



where, with now h{m) = g 



dt dm 
ym just involving the c-number m, 
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^ ^ dm 
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1 — m coth 



h{m) 



(17) 



(18) 



yh(m) ( h{m) 
w (m) = ^ ' I coth - ^ ' 

h 



(19) 



This type of equation has been extensively studied [10, 11]. We analyse below its solution in the present context. 
The term in 1/A^ of (O is negligible for smooth probabilities, but not for sharply peaked functions P{m), which 
occur at least in the initial state and in the final equilibrium ferromagnetic states. In fact, this term is dominant in the 
vicinity of the points where m = tanhh/T for which v(m) vanishes. In particular, as implied by the inequality il6i . 
P {m, t) tends for large times to an equilibrium shape characterized by the vanishing of the square bracket in illi . This 
condition entails 

N / m — mp ^ ^ 
2 



P{m, 



exp 



(20) 



where mp and 5p are given by 
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(21) 



(22) 



As readily checked, we recover dynamically the expected canonical distribution, including corrections for finite but 
large A^. However, for g ^ Jm^, eq. ( 12 H has two solutions such that ( I22t is positive, with nearly opposite values imp. 
Thus the asymptotic limit of P {m, t) for large t is not a priori known: it is a linear combination of the two ferromagnetic 
distributions MQ\ peaked around +mp and — mp. The present dynamical approach is necessary to explain how the 



evolution towards equilibrium produces the breaking of invariance. In fact, the weights ^+ and of the two terms 
+mp and —nip can only be found by solving the dynamical equation (I17> . which will allow us to express them in 
terms of the initial condition P (m,0). In the measurement problem, it is essential that in the If sector, for g > 0, the 
probability P (m,f ) concentrates for large times only around the positive ferromagnetic value +nip; a finite weight 
for the peak —nip would mean a wrong indication of the apparatus for the system in the state 5, = + 1 . Moreover, we 
need the lifetime of the initial paramagnetic state to be sufficiently large so that the interaction between S and A can 
be turned on while M still lies in the metastable paramagnetic state. In the theory of phase transitions, we wish to find 
through which intermediate shapes P {m, t) passes when the magnet M relaxes from the initial state to either one of the 
ferromagnetic equilibrium states, and to determine the probabilities of both occurrences. 

Our equation Mil , which governs the time-dependence of the order parameter, treated as a random variable, has the 
same form as the Fokker-Planck equation for a Brownian particle in one dimension submitted to an external influence. 
Its first term represents a deterministic drift, its second term a diffusion process which tends to widen the distribution 
P{m,t). Let us first drop this diffusion term. Eq. illl has then elementary solutions of the form 5 (m — m (/x,f )), where 
m (/i , f ) is the trajectory of a particle with initial position m — jj. and (position-dependent) velocity dm /dt~v (m) . (For 
very short times, which are not relevant below, we find from (I12> . ( I15> that v also depends on time, as v ~ — ' yP^fm, 
the sign of which implies that the fixed point m = is initially stable.) The functions m (/i , f ) is obtained by inverting 
the equation 

_ dm' _ h dm' 

~ J^i v{m')~YJ^ h{m')[l-m'cothh{m')/T] ' 

This motion of m has three fixed points, the zeroes of v (m). The closest to the origin, which for g <^ Jm^ lies at 

'«p = -j3y, (24) 

is repulsive: if ji differs slightly from mp, the point m{ii,t) moves astray. The other two fixed points are attractors 
associated with the ferromagnetic phases near +mp and —nip. Thus, if the diffusion term is discarded, the initial 
probability P (m,0) is split into two parts ni > nip and ni < nip, which will eventually result into sharp peaks located 
at +nip and —nip, respectively. In the evolution of P{ni,t), the weight P {ni, t) dm is conserved along the motion ( I23> . 
Hence, denoting by jx {m,t) the inverse mapping of m {lJ.,t), the solution of eq. (I17t without the diffusion term is [11] 

P(m,t)=P[^{m,t),0]'^^t^pp^, (25) 
V (ni) 

since dm/ dpi — v{m) /v(/x) forgiven t. We denoted v{ii{m,t)) as v[^(m,f)] and (m,f) ,0) as P[pL (m,f) ,0]. 

However, as we already noted, diffusion is essential in the first stage of the motion, when P{m,0) is concentrated 
around small values of m where v (m) is small. In this region, it is possible to solve the full equation (I17> by expanding 
V (m) sa (7/^) [g+ {J — T)ni] and w [m) k, jT /h'm powers of m, leading to 

h dp d T d^P 

7 at dm N dni-^ 

We can then find the exphcit solution of J23> . and invert the mapping m{pL,t), 

m{ii,t)^^e"' + j^(e"^-\), ^ (,n,f) = - ^ (l , (27) 



where we introduced the time-scale 
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e = , ^, . (28) 



rij-T) 

The solution of j26> , found by means of a Fourier transform on m, takes after some calculations the form 



where 



PM = / .^.-^V2Cp[^(^^,) + ^^o]^:[^ , (29) 



C^[l-e-^'/'^)jL-. (30) 



The result ( I29> encompasses the drift induced by the first term of i26\ and the diffusion induced by its second term. It 
turns out that the cumulated effect of diffusion is equivalent to a blurring of the initial value from which m is issued in 
the deterministic motion i27\ . over a width Vc which increases with t. We note however that, for e'^^ ^ 1 the effect 
of diffusion remains constant, as C tends then to Ccxd ~ T / [J —T). 

Since, according to \21\ . 6 is also the characteristic time over which m{jj.,t) exponentially diverges from its initial 
value we can match the effect ( I29l l of diffusion with the effect ( 125 1 of drift, even in regions where m is no longer 
small, by writing the solution of (I17> in the form ( I29> where {m,t) is now obtained from (12 3> rather than from ( I27t . 
This result holds, up to the very large times when the expression ( I29l l gets concentrated near ±mp. In the later stage of 
the evolution, that we do not need to consider, the diffusion term becomes again effective, and it determines the shape 
and width of the final probability according to ( I20> -(I22>. In order to discuss the behavior of ( I29t we note that when 
7 — r ^ 7, mp ~ 3 (/ — T) /J, one has mp ^ —g/ {J —T) nip, while v (m) has the form 



V (m) 



(in — mp) OTp — ) 



e 



(31) 



More generally, we note that \2>l\ has for g <^ Jin^ but arbitrary T < J the same features as dlSt : same zeroes, same 
behavior for small m. For this qualitatively good model of v (m) we can integrate explicitly (I23> for \mp\ ^nip as 



ni — nip I nip — jl^ 
jl—nip]] nip — m?- 



(32) 



or equivalently 



m{n,t) 



[pLe'l''~inp{e'l''-\)] 



jX {in,t) = nip + 



We take as initial state the narrow distribution 

P{in,0) = 



nip 



^nil+{^-nipf{e^'/S-l) ' 
(m — nip) e^'l^nip 



72-m2(l-e-2'/6 



/ N 



exp 



/ m — niQ 
2 I 5q 



(33) 



(34) 



(35) 



where otq characterizes a possible deviation from the paramagnetic state (for which otq — 0). The width Ss, determined 
by the initial equilibrium of M at a temperature To higher than J, is given by 5q =Tq/ {Tq — J). When quenching from 
Tq = oo, as is done in the figures, one has niQ —Q and 5o = 1 • Use of ( I32t in ( I29t and integration yield 
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27t{C + 8l 
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(36) 



where C is given by ( I30> . jJ. is the function of ni and f given by ( I34> . and dpi /dm = v (/x) /v{m). 

The expression ( I36> holds at all times, except in the very last stage of the evolution, when f (m,f), already 
concentrated near +mp and —nip, switches into the shape ( I20> . We can use it, however, to evaluate the probability 
or that the initial state (I35> ends up at one or the other ferromagnetic states +nip or —nip. Indeed, in the 
mapping m(/i,f), the points pL associated with the vicinity of the final point ni = +nip are those for which pL > nip. 
Taking thus ji instead of in as an integration variable in (I36> and letting f ^ oo, we find 



^_ = ierfc 



2 5 



(37) 



where 



-mp + niQ - 



J-T 



niQ 



(38) 



measures the bias due to the external field g and to the shift niQ in the initial distribution, where 



Co. + 5^ = j-j: + 5o' (39) 



accounts for the modification of the initial width 5o / ^/N due to diffusion, and where 



erfc(.^:) = I e''^ dt . (40) 
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UNIQUE FINAL STATE 

We can therefore distinguish two regimes which display several qualitatively different features. In the first case, the 
argument of erfc in ( I37> is large, so that the magnet M nearly certainly reaches the ferromagnetic state +m-p: the 
bifurcation is ineffective. This occurs, for an unbiased initial state (mo — 0), when 

{J-T)5 

^»^-=^; (41) 



this also occurs for g — Q when the shift mo of the initial state with respect to paramagnetism satisfies 

5 

mo > -= . (42) 



In this situation, the duration of the relaxation process is of order Q —h/y{J—T). More precisely, if we consider that 
ferromagnetism is reached when the distribution P{m,t) is peaked near m = 0.95 mp, this time, the registration time 
for a measurement, is found from i32\ to be 

In— ^. (43) 



'^^ 7iJ~T) b 

Since the range f < Treg does not involve A^, it is seen from ( I34t that A' occurs in the expression i36\ of P {m,t) only 
as a factor in the exponent. Hence P {m,t) displays at all times a single peak, see Fig. 1, that is narrow as 1/ \/N and 
located at the point m where fi {m,t) equals mo. The position of this peak is thus given by ( I33> in which /i is replaced 
by mo. Its width is seen from ( I36> to equal 5 (f) /Vn, where 

5„)^V^i|i|^ (44, 

When m(mo,f) is still close to the origin, diffusion lets C increase up to Coo = T / (J — T). The widening of the 
distribution P{m,t) is then governed by the dispersion of the speed v(m) of the drift motion [10]: the head of the 
distribution, with larger m, proceeds faster than its tail. The maximum of 5 (t) is thus attained when the peak of 
P{m,t) is located at the value of m which lets v(m) be maximum, that is, m — reached for 

h mp 

f = — rln^^. (45) 

We have then 

Imp 5 

&,„ = 5^. (46) 

Later on, 5 (f ) decreases, again due to the gradient in the drift velocity, which is now negative. This effect is, in the 
last stage, counterbalanced by diffusion which prevents 5 {t) from decreasing below the equilibrium value 8p given by 
(|23. 
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FIGURE 1. Transition from the initial paramagnetic state to the ferromagnetic state with magnetization up, driven by the couphng 
g. We have solved eq. (15) in the short-memory approximation, valid for y <C 1, where K,{co) is replaced by (7), the Debye cutoff 
being irrelevant. The parameters were chosen as A' = 1000, T = 0.657, g = 0.057, small but satisfying the condition (41). The 
curves describe the probability distribution P{m,t) at the times (from left to right) f/6 = 0, 0.5, 1, 2.25, 3, 4 and 5; the latter curve 
is near the infinite time limit. As initial distribution we took (35) with mo = and ^ = 1, corresponding to the paramagnet Tq = oo. 
The center of the peak moves from m = to rrif = 0.89707 according to (23); its width first increases up to (46) then decreases 
down to (22). 



We now consider the second regime, for which the bifurcation is active in the sense that neither 3^+ nor vanish. 
Such a possibility for a macroscopic system to transit towards several (here two) possible final states is known as 
spinodal decomposition. This situation is a nuisance for the measurement problem, since it corresponds to a failure 
in the establishment of correlations between the indication of A and the residual state of S. For the phase transition 
problem, it corresponds to dynamics for which the external field g and the initial shift niQ are too small to fully 
determine the phase towards which the system A will relax. Since at large times P {m,t) must have two sharp peaks, 
we know that some samples will reach +mp, other ones — mp, but we are interested here in the history of these 
samples. A natural question is the following: do the samples that will end up at +mp behave as above, with a well- 
defined trajectory m {t) within fluctuations of order 1 / ^/N1 More precisely, how does the distribution P (w, f ) evolve 
before it displays its final two-peaked equilibrium shape? 

To answer this question, which arises when the conditions (I4H . J42> are violated, we assume that the bias due to 
both the external field and the initial condition is small for large A^ as 



where 5 is defined by ( I39t and where A is either finite or small. The probability distribution, given by i36\ . begins to 
widen under the effect of diffusion as 



where C increases according to ( l3Qt . Even in case the initial distribution i35\ extends only over values larger than 
mp, that is, for 5o -^ib — niQ — mp, ( I48> develops a tail which may extend to the region m < nip, and which will give 
rise to the non-vanishing probability i37\ . This initial widening is relayed by the widening due to the gradient of the 
drift velocities, which will have here dramatic consequences. Indeed, whereas in the first regime the whole relaxation 
process was achieved after a delay ( I43> independent of A^, the probability is now not yet concentrated near imp for 
times such that e'/^ is of order a/^V, see Fig. 2. In such a range of times, the distribution (I36> takes the form 



ACTIVE BIFURCATION 
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where we introduced Suzuki's scaling variable [11] 



The large limit is thus singular. During a long lapse of time, the probability distribution ( I49l l presents no narrow 
peaks in 1 / \/N, although it behaves originally as i35\ and ends up as a sum of two contributions ilOi around +mF and 
—mp. Instead, it is a smooth function, independent of and extending over a wide range of values of m. 

Let us illustrate this behavior by describing the evolution of ( I49> in the completely unbiased case g — niQ ~ X ~ 0, 
for which P (m,t) is symmetric (Fig. 2). The problem is formally similar to the so-called laser model [11] The initial 
peak j35> with mo = begins to widen and to reach a finite width. It progressively covers most of the interval —nip, 
+mp. In particular, at the time 

^ = 01nfe^/f^ (51) 

such that = I, the distribution 
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P{m,t)^J^ '"'^ exp(-^— ^ — j] (52) 
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is very flat: near m ~ 0, it behaves as P °^ exp (— 3OT^/4mp); the ratio P{m,t) /P{0,t) is still equal to 0.93 for 
m = O.Smp, to 0.84 for m = Q.6mp, to 0.65 for m = O.lmp. At later times, the origin becomes a minimum and two 
maxima appear, lying at 



2a2 / Ne-^'l^ml 
m = ±mp\^\ —^±mp\i\ — £- . (53) 

These maxima move away from the origin towards ±mp. Originally not pronounced, they get narrower and narrower. 
For sufficiently large times, when a becomes small, they turn into unsymmetrical peaks near mp (and —mp), having 
the shape 

, , , e^^l^dx 2 mp — m 
Pmt/m--— — , x=—j— , (54) 

with a maximum at x = |. Both peaks eventually reach the symmetric equilibrium shape \2Q\ when the diffusion term 
dominates again. We estimate the relaxation time by evaluating after which delay these two peaks arrive at ±0.95mF, 
which yields 

This quantity is much larger than the relaxation times which occur in the first regime, by a factor ln\/^V. The relaxation 
time becomes even infinite when N °°. 

For finite values of the bias, as measured by A defined by ( l39t . (l47t . the curve ( l49l l is no longer symmetiic, 
but the results are similar. For A > 0, the maximum initially at the origin moves as m ~ Xmp/a while widening. 
At the time J51> . there is still a single maximum, but the distribution has become very wide, as for A = 0, with 
larger values for m > than for m < 0. A second maximum appears later on, at a negative value mj of m given by 



— — Amp ymp — 2m2 / V6, and at a time such that a} — ^ {m^ — mj) (mp — 2m2) /mp. Here again the two maxima 
end up as narrow peaks near ±mp, but equilibrium is reached for +mp earlier than for — mp; the geometric mean of 
these relaxation times is given by \55l . The main difference with the symmetric case is the occurrence of unequal 
probabilities ^+ and for the final phases, given by ( I37> . which reads 

^± = ^erfc(TA) . (56) 



These values, determined by the dynamics, are expressed in terms of the various parameters of the problem and of the 
initial conditions through ( l39t , WT\ . 
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FIGURE 2. Transition from the initial paramagnetic state to the mixture of ferromagnetic states for the same situation as in Fig. 
1, A' = 1000, T/J = 0.65 but g = 0. The times are: t/O = 0, 0.5, 1, 2.25, 3, 4, 5 and 10 (from top to bottom at m = 0; the last curve 
is very near the infinite-time limit). The distribution, initially peaked as l/^/N, becomes flat in Suzuki's scaling regime (49)-(50) 
(with A = here). For large t, two peaks build up near nif = 0.87206 and — mf, with widths governed by the diffusion term ~ 1 /N. 

BURIDAN'S ASS EFFECT 

Many among the properties exhibited above have been encountered in different contexts long ago [10, 11]. Indeed, the 
equation (I17> which governs the dynamics of the present phase transition has a very general form, describing directed 
Brownian motion with an unstable fixed point. We have stressed the existence for large of two contrasted regimes 
for the behavior of the probability P{m,t), depending on the initial conditions and on the field. We propose to term the 
second, anomalous situation "Buridan's ass effect". According to an argument attributed to Jean Buridan, a dialectician 
of the first half of the XlVth century, an ass placed equidistantly from two identical bales of hay will stay there and 
starve to death because it has no causal reason to choose one or the other [12]. 

Here, likewise, in a symmetric situation where the initial bifurcation may lead the magnetization m to +mp or —nip 
with equal probabilities ^+ — S^^ = j, an infinite time ( I55> is required in the macroscopic limit N ^ °° for the 
system to relax. We find more generally the same huge duration of the relaxation, from the paramagnetic state to either 
one of the ferromagnetic states, when there is a bias which is sufficiently small so that neither one of the probabilities 
^+ and for reaching +mp and —nip vanishes. In such a situation, if is large but finite, "the ass" ni will have 
reached one or the other bale after the delay j55L and we can predict with which probability: if we perform the same 
experiment with many asses, we know which proportion will have attained either side. However, for times t < 6 the 
asses have practically not yet moved, whereas for times such that e'/^ is of order \/N, we can make no prediction: the 
asses are scattered on a wide range between the two bales. This dispersion of the values of m during a long delay is a 
characteristic feature of the phenomenon. 

In contrast, the situation is trivial for a larger bias. This may take place either under the condition ( 14 H when the 
field g is sufficiently large (a wind which pushes the asses), or under the condition (I42> if the initial state has a residual 
magnetization (the asses are significantly closer to one of the bales). The process then takes a finite time j43> when 
N and we can predict with little error (in 1 / ^/N) where are the asses at each time. 

The giant fluctuations which occur in Buridan's ass effect may be regarded as a dynamical counterpart of those 
which occur at equilibrium near a critical point. In both cases, the order parameter, although macroscopic, is not a 
well-defined quantity even in the large limit, so that its behavior must be described by means of statistical mechanics 
rather than standard thermodynamics. In the present situation, no temperature can be defined for the system M, during 
the irreversible process which leads it from a high-temperature equilibrium paramagnetic state (I35> to a ferromagnetic 
equilibrium (I20> at the temperature T of the bath, fixed below the critical temperature J. However we may argue 
that in some sense M crosses its critical point during the dynamical process, in case it has some probability to reach 
both phases +mp and —nip in the end. The well-known critical fluctuations and critical slowing down thus manifest 
themselves here by the large uncertainty about m shown by P{m,t) during the long delay i55\ . Giant fluctuations in 
the dynamics are also well-known for nucleation processes. However in the present long-range model there is no 
space structure, so that we have only to deal with the statistics at each time of the single order parameter m. 



CONCLUSION: BURIDAN'S ASS AND QUANTUM MEASUREMENT 



As regards the use of A = M + B as an apparatus to measure the spin S, two conditions should be satisfied. On the 
one hand, the initial metastable paramagnetic state of M, prepared by setting its temperature at Tq > J, should have a 
long lifetime i55\ . so that the interaction Q between S and A can be turned on before M has begun to relax. This is 
achieved provided there is neither a significant external field go nor a lack of symmetry mo before the measurement, a 
condition expressed more precisely by A ^ 1 in (I47> . that is, using 5q =To/ {To — J), 

SO , V ^ f , ^0 \ 

' mo < 77 1^ — ^ + ; • (57) 



J~T J N\J-T To-J 

On the other hand, the coupling g between S and A should be sufficiently large so as to ensure faithfulness of the 
registration, +mp in the sector ft, —mp in the sector j j. This is achieved provided Buridan's ass effect does not take 
place during the measurement process, that is, provided 

■»i('^ + ^V (58, 



J-TJ N\J-T To-J 

Then the registration time j43> is finite for large A^. 

We have chosen to study the model (0} because its exact solution at equilibrium can be obtained for large through 
a static mean-field approach; we had the prejudice that its dynamics might also be solved exactly through a time- 
dependent mean-field approach. It turns out that, although an exact solution in the considered limiting case is available, 
it cannot reduce to mean-field when Buridan's ass effect is present. A time-dependent mean-field approach relies on 
the existence of trajectories m{t) around which m has negligible statistical fluctuations. However, when the bias b 
defined by ( I47> is small, the dynamical process is governed by fluctuations which become macroscopic. Even if the 
initial magnetization is exactly defined {5o = 0) the diffusion term in illi produces a width j39> which eventually 
leads to a flat distribution, as illustrated by i52\ . It is only in the biased regime, which leads the magnetization to a 
single value +mp with nearly unit probability, that fluctuations play little role for large A^, and that the time-dependent 
mean-field equation i23\ is sufficient to describe the dynamics. When the bifurcation is active, the intuitive idea that 
the variable m, because it is macroscopic, should display relative fluctuations small as 1 /\/^V becomes wrong except 
near the initial time and near the final equilibrium states. 

The two regimes are characterized by the final situation: either non-vanishing probabilities £P+ and for the two 
alternatives +mp and —mp, or a single choice only. Everything takes place, like for Buridan's ass, as if the process 
were governed by final causes: the behavior is deterministic if the target is unique; it displays large uncertainties in 
the dynamics if hesitation may lead to one target or the other. In fact, the singularity of the problem arises from the 
structure of ( I17> . which contains a deterministic drift term involving the velocity v (m) and a diffusion term. Small 
as the latter term does not contribute significantly when the bias is sufficient to determine a single outcome 
(|A| ^ 1); it becomes essential when the bifurcation is active {X finite or smafl), because it governs the right-hand 
side of jl7l l near the fixed points of the drift motion, where v (m) vanishes. The long duration of the relaxation and the 
large uncertainties in its dynamics then reflect in a probabilistic language the slowness of the pure drift motion around 
the unstable fixed point, and the long and random delay needed to set m into motion. The very direction of the drift, 
which will be reflected in the sign of the final magnetization, is thus quite sensitive to the perturbation caused by the 
diffusion term, which altogether determines the long-time behavior. 
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